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(57) Abstract: The invention concerns a method for packing and filtering 
geophysical events read on multi-domain recordings, with distribution of 
said events into families each having a particular geophysical significance: 
equal offset or equal -angle incidence data cubes, elastic parameter cubes 
derived from a joint stratigraphic inversion and the like, in order to ex- 
tract therefrom data on the type of subsoil, using a multivariate statistical 
technique. The method essentially comprises forming by combination of 
geophysical variables, synthetic variables much fewer in number, which 
are obtained by constructing an orthogonal vector base in each of the sets 
of analysis constituted by the data of each of the families, wherefrom is de- 
rived the formation of an orthonomal vector base (new attributes) for filter- 
ing and describing said geophysical events. The invention is applicable to 
geological interpretation or to interpretation of an underground reservoir, 
from seismic measurements or diagraphic measurements, for example. 

(57) Abrege : M6thode pour compacter et filtrer des evenements gSophy- 
siques lus sur des enregistrements multidomaines, avec repartition de ces 
evenements en families avec chacune une signification physique particu- 
liere: cube de donnees iso-d6port ou iso-angle d' incidence, cubes de para- 
metres elastiques issus d'une inversion stratigraphique jointe, etc. dans le 
but d'en extraire des informations fur la nature du sous-sol, en exploitant 
une technique statistique multivariable. La m6thode comporte essentiellement la formation par combinaison des variables g6ophy- 
siques, de variables synthStiques en nombre tres infeneur, que Ton obtient par construction d'une base vectorielle orthogonale dans 
chacun des ensembles d'analyse constituSs par les donnees de chacune des families, d'ou decoule la formation d'une base vecto- 
rielle orthonormee d^crivant ces ensembles d'analyse, et 1' utilisation de cette base vectorielle orthonormee (nouveaux attributs) pour 
filtrer et decrire les dits 6v6nements geophysiques. Applications a 1' interpretation g£ologique ou a l'interpr6tation d'un reservoir 
souterrain, a partir de mesures sismiques ou de mesures diagraphiques, par exemple. 
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5 METHODE POUR COMPACTER ET FELTRER DES EVENEMENTS 
GEOPHYSIQUES DANS LE BUT D'EN EXTRAIRE DES INFORMATIONS SUR 
LA NATURE DU SOUS-SOL 

La prSsente invention concerne une methode pour compacter et filtxer des donn6es 
10 g6ophysiques dans le but d'en extraire au mieux des informations sur la nature du sous-sol. 

Cette m6thodologie s' applique k la compaction et au filtrage des jeux de donn6es 
geophysiques multidomaines, i.e. soit des donnees sismiques avant sommation issues 
d'une campagne sismique 3D regroup6es en un certain nombre de cubes de donn6es 
correspondant chacun par exemple k un meme deport (cubes de traces sismiques dits iso- 

15 offsets) ou k un meme angle ^incidence (dits cubes iso-angles d'incidence), ou k des cubes 
d'attributs sismiques obtenus par traitement des amplitudes sismiques comme, par 
exemple, des cubes de propri6tes 61astiques issus d'une inversion stratigraphique jointe de 
donnees sismiques avant sommation, ces cubes etant obtenus par la mise en ceuvre de la 
methode decrite par exemple dans le brevet FR 2 800 473 (EP 1 096 270) du demandeur. II 

20 s'agit alors de mettre en forme Finformation sismique avant sommation contenue dans des 
portions de traces sismiques form6es par exploration d ! une zone souterraine, h. l'aide de 
nouveaux attributs repr6sentatifs qui sont ensuite exploits pour 1' interpretation g6ologique 
ou Tinterpr6tation d'un reservoir du sous-sol. Elle peut aussi s'appliquer k la compaction et 
au filtrage de donn6es sismiques 3D avant ou apr6s sommation de campagnes r6p6t6es 

25 dans le temps (sismique 4D), ceci afin d*interpr6ter les variations de contenu en fluides 
d'un reservoir p&rolier en cours de production. Cette m6thode peut aussi s ! appliquer k la 
compaction de donnees diagraphiques de plusieurs natures physiques, en formant k partir 
des mesures de depart de nouveaux attributs diagraphiques, et ce afin d'ameliorer 
Interpretation des donnees diagraphiques en termes de caract6ristiques lithologiques et 

30 p^trophysiques du reservoir. 
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Etat de la technique 

Les mesures geophysiques sont classiquement utilisees pour fournir des informations 
compl6mentaires, par rapport aux observations directes dans les puits, sur les .variations des 
propriet6s g6ologjques des formations du sous-sol. L' interpretation sismique, par exemple, 
5 est bas6e sur r analyse d'attributs afin d'extraire ces proprietes gdologiques k partir de 
l'information sismique, et ce dans une zone d'int^ret d'une formation souterraine (zone 
r6servoir par exemple). Un autre exemple est Interpretation diagraphique ou les diverses 
mesures sont analysees simultanement pour en deduire les caracteristiques p6trophysiques 
de la formation souterraine etudi6e. 

10 Une approche standard pour T analyse d'attributs sismiques ou diagraphiques 

implique r utilisation de m6thodes statistiques de reconnaissance de forme et d* estimation 
afin de calibrer les attributs pertinents avec rinformation disponible obtenue au droit de 
diff6rents puits for6s. 

Par exemple, par le brevet FR -A 2 768 818 (US 6 051 651 du demandeur, on connait 
15 une m&hode permettant de d6tecter les classes potentielles dans une population 
d'evenements li6s aux propri6t6s physiques d*un milieu complexe tel que le sous-sol, 
rep6r6s k partir de donnSes obtenues par exploration du sous-sol ou par des mesures in situ, 
(6venements d f ordre g6ologique, g6ophysique, ev6nements li6s k des donn6es de 
production de fluides etc.), ces 6venements etant reper6s par des points dans un espace 
20 multivariable et definis par un vecteur d'attributs. 

Un autre exemple est donn6 par le brevet FR 01/05.675 du demandeur, qui d6crit une 
m6thode pour faciliter le suivi au cours du temps de revolution d'une zone souterraine par 
analyse compare d'un certain nombre n de jeux d'enregistrements sismiques obtenus 
successivement k Tissue de n campagnes sismiques 3D successives (sismique dite 4D), 
25 dans laquelle on utilise une technique de reconnaissance des formes appliqu6e k Tensemble 
des 6v6nements sismiques de plusieurs campagnes, considers et analyses simultanement. 

Citons aussi le brevet EP -A-671 017 (US 5 638 269) du demandeur, qui d6crit une 
m6thode permettant .de rapprocher des donn6es gdologiques obtenues par carottage ou 
diagraphies dans des puits et des donndes sismiques obtenues par une campagne 
30 d'exploration sismique qui s'appuie sur une technique de calibrage statistique avec un 
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etalonnage statistique obtenu en rapprochant les donn6es g6ologiques locales (mesurees 
dans des puits) avec des attributs sismiques lus sur des traces sismiques obtenues au 
voisinage imm6diat de chaque puits. 

Les attributs sismiques qui font l'objet de ce type ^interpretation sont calcules dans 
5 le domaine apres sommation, ou dans le domaine avant sommation. 

Dans le domaine apres sommation, les attributs classiques sont calcules a partir des 
amplitudes au niveau du reservoir, ou k partir de Timp6dance P estimee k l'aide d'une 
inversion stratigraphique. 

Dans le domaine avant sommation, le nombre d' attributs pertinents peut 
10 consid6rablement augmenter. Dans ce cas, les attributs sont soit les amplitudes a differents 
offsets (d6ports) ou angles d'incidence, soit des parametres issus d'une inversion 
. stratigraphique jointe avant sommation (Impedances P et S, masse sp6cifique, produit des 
parametres de Lam6 avec la masse sp6cifique, etc.). 

En ce qui concerne rinterpr6tation des mesures diagraphiques, la situation est 
15 analogue, un tr&s grand nombre de mesures 6tant disponibles d'oii en decoulent un grand 
nombre d'attributs permettant de caracteriser les diagraphies, ces attributs 6tant souvent 
organis6s en families k m§me signification physique (par exemple, famille des attributs de 
r6sistivit6, famille des attributs associ6s k la radioactivite naturelle ou induite des 
formations, etc.). 

20 Une m6thode dite danalyse en composantes principales (ACP) bien connue des gens 

de Tart, peut etre appliqu6e k ces attributs afin d* analyser les relations existantes entre eux 
et de r6duire le nombre d' attributs significatifs. Les composantes principales extraites de 
r ACP sont les nouveaux attributs : ils d6finissent une base vectorielle orthogonale ou non 
(rotation) et sont des combinaisons lin6aires de variables initiales. lis peuvent Stre utilisees 

25 dans des algorithmes statistiques de reconnaissance de forme et correspondent k un filtrage 
multivari6 de Tinformation sismique ou diagraphique initiale. Differentes applications de 
T ACP sont d6crites, par exemple, dans : 

- Dumay, J., Fournier, F., 1988, "Multivariate statistical analyses applied to seismic 
facies recognition", Geophysics, 53, 1151-1159 ; 
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- Hagen, D., C, 1982, "The application of principal components analysis to seismic data 
sets" : Geoexpl., 20, 93-111 ; ou encore 

- dans le brevet FR-A- 2 772 138 (US 6 345 108) du demandeur. 

Mais TACP ne prend pas en compte la notion de groupe d'attributs et, par 
5 consequent, il est souvent difficile de donner une interpretation physique aux composantes 
principales, ou bien de les relier clairement aux attributs initiaux, surtout si ceux-ci sont 
organis6s en groupes de meme signification physique. 

U analyse dite canonique est aussi une methode statistique bien connue des gens de 
Tart, qui permet d'6tablir les relations pouvant exister entre deux groupes de variables afin 
10 de savoir si ces deux groupes decrivent les mSmes propri6t6s. Cette methode est decrite par 
exemple par : 

- Hotelling, EL, 1936, "Relations between two sets of variables", Biometrika, 28, 321- 
377. 

Un exemple d' application de T analyse canonique est d6crite par exemple par : 

15 - Fournier, R, and Derain, J.F., 1995, "A statistical methodology for deriving reservoir 
properties from seismic data", Geophysics, 60, 1437-1450. 

Cette methode se limite k F6tude de deux groupes de variables et done ne s' applique 
pas de fagon generate k l'information sismique multicubes ou k 1'information diagraphique 
multidomaines. De plus, elle definit des variables synthetiques dans chacun des sous- 
20 espaces associ6s aux deux groupes, et non une seule base vectorielle permettant de 
compacter Tensemble des variables initiales et de d6crire les relations entre les deux 
groupes. 

Plusieurs methodes permettant de g6n6raliser 1' analyse dite canonique en 6tendant 
T analyse k plus de deux groupes de variables existent. Diff6rents aspects de Tart anterieur 
25 dans le domaine considere sont ddcrits, par exemple, dans les publications suivantes : 

- Horst, P., 1961, "Relations among m sets of measures", Psychometrika, 26, n°2, 129- 
149 ; 
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- Caroll, LD., 1968, "A generalisation of canonical correlation analysis to three or more 
sets of variables", Proc. 76 th Conv. Amer. Psych. Ass ; 

- Kettenring, J.R., .1971, "Canonical analysis of several sets of variables", 
58, 3, 433-450 ; ou 

5 - Saporta, G., 1990, "Probability, analyse des donnees et statistiques", Technip, Paris. 

Ces methodes presentent certaines limites. En particulier, elles ne permettent pas de 
decrire chaque groupe separ6ment, les variables synth6tiques ne peuvent etre relives 
- facilement aux differents groupes, et , par consequent, etre physiquement interpr6t6es. De 
plus, seule la part d'inertie globale qu'elles reprSsentent peut etre connue, et non la part 
10 d'inertie qu'elles repr^sentent pour chaque groupe. H n'est par cons6quent pas possible de 
filtrer chaque groupe ind6pendamment. 

Une technique d'analyse 6galement connue dite analyse en composantes principales 
g6n6ralis6e (ACPG), pennet de comparer diff6rents groupes de variables tout en d6crivant 
chacun d'eux. On la trouve mise en oeuvre sur des donnees k caractdre 6conomique par 
15 exemple dans : 

- Casin, Ph, 2001, "A generalization of principal component analysis to K sets of 
variables", Computational Statistics & Data Analysis, 35, 417-428. 

La methode selon Pinvention 

20 La methode selon Pinvention propose une approche oh en exploitant une technique 

statistique multi variable appliquee & plusieurs jeux d'enregistrements sismiques ou 
diagraphiques disponibles, obtenus par exploration d'une zone souterraine, on compacte et 
on filtre les 6v6nements g^ophysiques lisibles sur ces jeux de mesures. La m6thode 
comporte : 

25 - la r6partition de ces 6v6nements g6ophysiques dans un certain nombre p de families 
d f attributs g6ophysiques ayant chacune une signification physique particuliere ; 



la formation par combinaison des dits attributs g6ophysiques, de variables synth6tiques 
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en nombre tres inftrieur au nombre des dits attributs, que Ton obtient par construction 
d'une base vectorielle orthogonale dans chacun des p ensembles d'analyse constitues 
par les donntes de chacune des dites families, d'oii decoule la formation d'une base 
vectorielle orthonormee dtcrivant les p ensembles d' analyse ; 

5 - Tutilisation de cette base vectorielle orthonormte pour filtrer et decrire les dits 
tvtnements gtophysiques. 

La formation des variables synthttiques est obtenue par exemple en combinant une 
technique d' analyse en composantes principales ou ses techniques dtrivtes bashes sur la 
diagonalisation d'une norme (recherche de vecteurs propres) et une technique de 
10 regression simple. 

Suivant un mode de mise en oeuvre, on peut soumettre les variables synthttiques a 
une rotation pour amtliorer les correlations entre variables initiates, c'est k dire les attributs 
geophysiques, et variables synthttiques. 

Les variables synthetiques peuvent etre form6es k partir de families comportant par 
15 exemple un ou plusieurs cubes de traces sismiques obtenues k Tissue d ! une campagne 
sismique 3D, chacun correspondant k un meme d6port ou offset. 

On peut egalement former des variables synthetiques k partir de families constitutes 
de cubes de traces sismiques obtenus k Tissue d'une campagne sismique 3D, chacun 
correspondant k un meme angle d'incidence. 

20 On peut encore former des variables synth&iques k partir de families constitutes de 

cubes d' attributs sismiques founts par traitement des amplitudes de traces sismiques 
obtenues k Tissue d'une campagne sismique 3D. 

On peut encore former des variables synthttiques a partir de families constitutes de 
cubes de traces sismiques obtenues k Tissue de campagnes sismiques successives de la 
25 zone. 

On peut encore former des variables synthttiques a partir de families constitutes de 
cubes d' attributs sismiques formts par traitement des amplitudes de traces sismiques 
obtenues k Tissue de campagnes sismiques successives de la zone. 
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On peut encore former des variables synthetiques a partir de families constitutes de 
mesures d ! attributs diagraphiques de nature physique differente, comme la famille des 
mesures de resistivit6, la famille des mesures liees k la radioactivity des terrains, la famille 
des mesures li6es k la porosit6 des terrains, etc. 

5 La m6thode est particulierement avantageuse en ce qu'elle permet, appliquee k de 

Tinformation sismique ou diagraphique: 

- de d6finir de nouveaux attributs g6ophysiques, combinaisons linSaires des attributs 
geophysiques initiaux, resumant Tinformation sismique ou diagraphique initiale et 
permettant de fiQtrer cette information, tout en respectant la structure multidomaines 

10 des donn6es 

- de montrer si toutes les families d' attributs geophysiques initiaux, ou seulement 
certaines d'entre elles, apportent une information redondante, et de mesurer la part de 
cette information redondante 

- d'apporter une interpretation physique aux nouveaux attributs g6ophysiques. 

15 Les nouveaux attributs g6ophysiques ainsi form6s peuvent ensuite Stre analyses afin 

d ! extraire les propri6t6s g6ologiques de Tinformation sismique ou diagraphique. Deux 
classes de m&hodes peuvent classiquement Stre utilises k cette fin : une premiere qui 
utilise des techniques d'estimation statistique et fournit une interpretation quantitative en 
termes de proprittfis gtologiques, et une deuxi&me qui utilise des techniques de 

20 reconnaissance de forme pour analyser les faci&s sismiques ou les faciSs diagraphiques 
(electrofacies) et fournir une interpretation qualitative. 

Des exemples de ces deux types de m6thode sont mis en ceuvre respectivement dans 
les brevets EP-A-671 017 (US 5 638 269) et FR-A-2 768 818 (US 6051651) du 
demandeur, d6jk citds. 

25 

Presentation succincte des figures 



les figures 1A k 1C montrent T organisation des enregistrements sismiques en p 
groupes : chaque groupe rassemble dans un cube T ensemble des portions de traces 
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correspondant h un crit&re physique identique ; 

- les figures 2A a 2C montrent respectivement les cartes de trois premiers nouveaux 
attributs sismiques Z 0 *, et Z (3) ; 

- les figures 3A k 3D montrent respectivement le premier plan temporel vrai (3A) de la 
5 fenetre d'analyse et le meme plan reconstituS/filtre avec 1 (3B), 3 (3C)ou 8 (3D) des 

nouveaux attributs sismiques pour le cube d' impedance P ; 

- les figures 4A k 4D montrent respectivement le premier plan temporel vrai (4 A) de la 
fenetre d' analyse et le mSme plan reconstitue/filtr6 avec 1 (4B), 3 (4C)ou 8 (4D) des 
nouveaux attributs sismiques pour le cube d' impedance S ; 

10 - les figures 5 A h 5D montrent respectivement le premier plan temporel vrai (5 A) et le 
meme plan reconstitu6/filtre avec 1 (5B), 3 (5C)ou 8 (5D) des nouveaux attributs 
sismiques pour le cube de masse sp6cifique p ; 

- les figures 6 montrent respectivement les sections verticales, passant par le puits W2, 
vraie (6 A), reconstituee/filtree avec les 8 premiers nouveaux attributs (6B), ainsi que 

15 les differences (r6sidus) (6C) entre ces deux sections pour le cube d'impedance P ; 

- -les figures 7 montrent respectivement les sections verticales, passant par le puits W2, 
vraie (7 A), reconstitu6e/filtree avec les 8 premiers nouveaux attributs (7B), ainsi que 
les differences (r6sidus) (7C) entre ces deux sections pour le cube d' impedance S ; 

- la figure 8 montre la localisation des traces d'apprentissage utilis6es pour calibrer la 
20 fonction de classement dans le cadre de la reconnaissance des formes supervisee des 8 

premiers nouveaux attributs sismiques ; 

- les figures 9A, 9B montrent respectivement les attributions des 6v6nements sismiques 
et les attributions les plus fiables de ces memes ev6nements dans le cadre de la 
reconnaissance des foimes supervis6e des 8 premiers nouveaux attributs sismiques ; et 

25 - les figures 10A, 10B represented respectivement les attributions des 6v6nements 
sismiques et les attributions de ces mSmes 6v6nements sismiques dans le cadre de la 
reconnaissance des formes non supervisee des 8 premiers nouveaux attributs sismiques. 
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Description detaillee 

On dispose de plusieurs jeux d'enregistrements geophysiques obtenus par exploration 
sismique ou reconnaissance diagraphique d'une zone souterraine etudi6e. Les 6venements 
geophysiques que Ton cherche k compacter/filtrer sont caract6rises par- des attributs 

5 g6ophysiques souvent organises en groupes ou families a signification physique 
particuli£re (figures 1A k 1C). Ces attributs peuvent etre la succession des valeurs des 
amplitudes de portions de traces sismiques ou de paramdtres 61astiques issus d'une 
inversion 61astique jointe avant sommation, obtenue par exemple a l'aide de la m6thode 
dScrite dans le brevet FR 2 800 473 (EP 1 096 270)d6jk cite le long de la fenetre 

10 temporelle 6tudi6e. Ces attributs peuvent etre, pour une profondeur donnee, les mesures 
d f un meme ensemble d'enregistrements diagraphiques r6alis6s le long d'un puits, comme 
par exemple les valeurs des diff6rents outils de r6sistivit6, les valeurs des differents outils 
mesurant la porosit6, etc. 

1 - Notations 

15 Soit X une matrice compos6e de p tableaux juxtaposes Xi,. . . Xj . . . X p : 

X = {Xj ... Xj ... X p } , Xj = [Xj j ... Xj k ... Xj m ^] 

oil Xj est un tableau compos6 de mj variables Xj, k de n individus (au sens statistique du 
terme). La matrice X rassemble tous les 6v6nements geophysiques disponibles. Ces 
6v6nements sont tout d'-abord organises en p groupes ou families Xj (figure 1A). Chaque 
famille Xj correspond a une propriety physique particuliere. Pour Tanalyse de la sismique 
20 3D, cette famille se pr6sente sous la forme d ! un cube de valeurs correspondant a cette 
propri&6, rep6r6es par leur position spatiale (ligne, cdp, temps/profondeur), (figures IB et 
1C). Ces cubes peuvent etre form6s par exemple par un ensemble des valeurs d'amplitude 
de n traces ou portions de traces associees k un deport commun ou k un angle commun ou 
encore la distribution spatiale des valeurs d'un attribut commun d6riv6 de ces amplitudes 
25 telle que, par exemple, l'imp6dance acoustique en ondes P ou S, la masse specifique ou des 
param&tres de Lam6, dans une fenetre temporelle choisie de taille mj. Dans cet exemple, 
chaque variable Xjjc rassemble les valeurs de V attribut sismique consid£r6 (amplitude 
sismique ou propri&6 61astique) pour un 6chantillon temporel k rint6rieur de la fenetre 
d' analyse. Dans cet exemple, un cube comprend autant d'individus que de traces 
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disponibles. Pour Tanalyse des enregistiements diagraphiques, cette famille se presente 
sous la forme d'une suite de series de valeurs, rep6r6es par leur position le long du puits 
(temps ou profondeur). Cette suite peut etre formee, par exemple, par l'ensemble des 
mesures diagraphiques des outils de meme nature physique, comme les diff6rents outils 

5 mesurant la porosite de la formation, ou les outils mesurant differents types de resistivite, 
sur un ensemble de n cotes profondeur ou temps le long d*un ou de plusieurs puits au 
niveau de la formation souteiraine etudi^e. Dans pet exemple, chaque variable 
rassemble les valeurs de Tattribut diagraphique considere (resistivite, porosite, 
radioactivite, ...) pour une cote profondeur ou temps k Tinterieur de la fenetre d* analyse. 

10 Dans cet exemple, une suite comprend autant d'individus que de cotes profondeur ou 
temps disponibles le long du ou des puits analyses. 

Les mj variables Xj, k d'une table Xj definissent un .espace denot6 Sj et sont suppos6es 
centrees. Soit D n la matrice diagonale des poids des individus. 

Soit 11^ un op6rateur de projection orthogonale sur Tespace Sj defini par : 

is n j = x J (x t j D n x j )- 1 x' j D n 



2 - Calcul des nouveaux attributs 

Le but de 1* Analyse en Composantes Principales Gen6ralisee (ACPG) est k la fois de 
dgcrire les proximit6s entre les p families de variables. et de d6crire separ6ment chaque 

20 groupe gr^ce a des variables synth&iques qui sont des combinaisons lin6aires des variables 
initiales. Alors, ce qui est commun & tous les groupes ou families, ou & certains d' entre eux 
settlement, et par consequent ce qui est different, peut Stre mis en evidence. L'ACPG 
permet de r6duire le nombre de variables utilis6es pour decrire un jeu de donnees et permet 
un filtrage multivari6 tout en conservant la structure multidomaines des donnees initiales. 

25 Afin d'atteindre ce but, la premiere variable synth6tique calculee, que Ton nomme Z (1) , 
calcul6e doit Stre telle que : 

1. les p projections orthogonales Zj (1) de Z (1) sur les espaces Sj doivent etre le plus 
proche les unes des autres ; 
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2. Dans chaque groupe j la dispersion des n individus projetes sur Zj (1) est 
maximum. 



1. max £r 2 (ZP\Z«) 

2. max 2Cov 2 (Zf>^, v ) 

Ceci se traduit math6matiquement par les conditions suivantes : 

oil R 2 est le coefficient de correlation multiple et Zj (1> est la projection de Zf 1) sur Tespace 
5 Sjnormeekl. 

Ensuite, on cherche une nouvelle variable synthetique Z* 2) selon les memes criteres. 
On ajoute une condition d'orthogonalite entre les projections les variables Z (k) du meme 
groupe. Le processus s'arrSte quand une base orthogonale dans chaque espace Sj est 
10 trouv6e. Le nombre de variables Z (k) n'excede pas la dimension la plus grande des espaces 

Des m6thodes de rotation bien connues des gens de Tart, peuvent 8tre mises en ceuvre 
dans chacun des espaces ainsi trouv^s, afin d'optimiser les coefficients de correlations 
entre les nouvelles variables et les variables initiates. 

15 La premiere variable Z (1) est vecteur propre norme de la matrice de covariance 

XX l D n correspondant k la plus grande valeur propre. Z (1) est ensuite projete sur les p 
espaces Sj et, dans chaque espace, la regression des mj variables par Zj (1) est calcul£e. La 
matrice des p residus juxtaposes X/ 2) definit le nouvel espace d'etude. Ce nouvel espace 
correspond & Tespace initial duquel Tinfluence de Zj (1) a ete retire dans chaque groupe Sj. 

20 Les mj variables de la table Xj (2) definissent un sous espace Sj (2) de Sj. La seconde 
composante Z (2) est le vecteur propre norme de la matrice de covariance X^X^Dn des 
residus correspondant k la valeur propre la plus grande. Z (2) est projetee sur les p espaces 
Sj (2) . Et la regression des mj variables par Zj (2) est calcuiee dans chaque espace Sj (2) . La 
regression effectuee dans chaque espace Sj assure que la projection de Z (2) sur r espace Sj (2) 

25 est orthogonale h la projection de Z (1) sur Tespace Sj. 
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Le calcul des vecteurs propres est effectue par toute m6thode connue des gens de 

Tart. 

En pratique, la matrice XX l D n etant d'ordre n, n pouvant etre tr&s grand, a une etape 
k, on calcule le vecteur propre norme V 00 de la matrice X (k)t D n X (k) associe k la plus grande 
5 valeur propre. Le vecteur Zf® est ensuite simplement calculi par transformation lineaire 
Z^X^V^etnonn^ 

Chaque groupe peut etre reconstitu6 avec un petit nombre d'attributs. Les variables 
initiales du groupe j peuvent s'exprimer comme la somme d'une combinaison lineaire des 
variables Zj 00 et d'un residu. Si toutes les variables Zj 00 d'un groupe sont utilisees, le 
10 r6sidu est nul. Les coefficients de la combinaison lineaire ne sont autres que les 
coefficients calcul6s k Fissu de chaque regression dans les espaces Sj. 

La m6thode selon invention a pour propriete, entre autres, que les variables Zf® sont 
orfhogonales deux k deux. 

En vue de 1' interpretation des variables synthetiques calcuiees, on calcule ensuite 
15 pour chaque groupe j : 

1. le carr6 de la correlation entre la variable Z (k) et la variable 2j W du groupe j, 
indiquant le degr6 de proximite entre ces deux variables; 

2, le pourcentage d'inertie du groupe j explique par la variable Zj W . 

Les correlations entre les variables Z (k) et les variables initiales Xjj sont egalement 
20 calcuiees pour mettre en evidence d'eventuelles relations entre les Xjjk . Les correlations 
entre les variables Zj W et Xj^ d'un meme groupe j, ou les correlations entre variables Zj (k) 
pour un meme k peuvent egalement §tre calcuiees. 

Les nouvelles variables synthetiques ainsi creees constituent les nouveaux attributs 
geophysiques recherches. 

25 3 - Exemple d 'application 

La methode selon Tinvention a 6te appliquee k trois cubes de propri6tes eiastiques 
(impedance P notee Ip, impedance S notee Is et masse specifique notee p ) obtenus apres 
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inversion elastique jointe avant sommation k Taide de la m6thode d6crite dans le brevet FR 
2 800 473 d6ja cit6 d'enregistrements sismiques regroupes en cubes & angles d'incidence 
commun. Chaque propri6t6 elastique engendre un groupe de variables ou d'attributs 
sismiques. La fenetre temporelle d'analyse se compose de 19 £chantillons de signaux 

5 pr61ev6s sur les enregistrements sismiques. Nous consid£rons comme variables les valeurs 
de la propriety 61astique pour un echantillon temporel k l'interieur de la fenetre d'analyse. 
Chaque groupe est alors compost de 19 variables compos6es d'autant d'individus que de 
traces disponibles. Comme il y a trois groupes de variables, nous disposons au total de 57 
variables. Les trois cubes de paramStres ont 6t6 pr6alablement horizontalisfis par rapport a 

10 une surface de reference. 

Les variables synthetiques Z 00 obtenues resument les informations respectives des 
trois propri6t6s 61astiques, tout en prenant en compte les relations entre elles. 

3.1 - Analyse des relations entre les impedances P, S et la masse specifique p 

Le tableau 1 ci-dessous donne la proportion de variance en pourcentage expliquee 
15 par les projections des 8 premieres variables synthetiques Z 0) & Z (8) dans les trois 
groupes. H montre qu'il suffit de 8 variables pour resumer 87,4% de la variance de 
Fimp6dance P, 94,4% de la variance de Fimp6dance S et 88,4% de la variance de la masse 
sp6cifique. 





Z 0) 




2» 


Z (4) 


Z (5) 






2» 


• 


Ip 


23,4 


24,5 


6,3 


13,1 


9,8 


3,1 


3,8 


3,4 


87,4 


Is 


26,6 


14,2 


15,2 


12,6 


5,7 


7,2 


5,8 


7,1 


94,4 


• 


16,9 


144 


13,6 


10,2 


8,0 


9,4 


10,1 


6,1 


88,4 



20 

Le tableau 2 ci-apres donne le carr6 de la correlation entre les variables synthetiques 
Z w et leurs projections sur les trois groupes. Cette valeur indique le degr£ de proximite 
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-entre la variable et le groupe. Plus la valeur est proche de 1, plus la variable est proche de 
sa projection et done du groupe. 









z cs) 


Z (4) 


Z (5) 


Z (6) 


z cd 


Z (8) 


Ip 


0,48 


0,79 


0,38 


0,72 


0,71 


0,24 


0,02 


0,30 


Is 


0,91 


0,43 


0,86 


0,53 


0,36 


0,59 


0,09 


0,73 


• 


0,78 


0,49 


0,77 


0,46 


0,65 


0,79 


0,96 


0,46 



5 Par exemple, ZP est tres proche de Is et proche de p , mais tres 61oigne de Ip. Alors 

que Z5 2) est proche de Ip, mais tres 61oign6 de Is et de p. Sur les 8 premi&res variables Zf®, 
trois sont proches de Ip seul, deux sont proches a la fois de Is et de p , une est proche de Is 
seul, et deux sont proches de p seul. Ainsi, il apparait que Ip et Is ne contiennent pas la 
meme information et qu'une part de T information contenue dans Is et dans p semble 

10 commune. 

3.2 - Filtrage des impedances P, S et de la masse specifique 

Les trois propriet6s 61astiques peuvent etre reconstituees avec un nombre restreint de 
nouveaux attributs sismiques. Le nombre de variables synth6tiques utilises pour la 
reconstitution peut etre different d'un groupe & T autre : il depend du pourcentage d'inertie 
15 que Ton d6sire reconstituer/filtrer. 

3.3 - Reconnaissance des formes supervisee et non supervisee des nouveaux 
attributs 

II est possible d' interpreter geologiquement les nouveaux attributs sismiques au 
travers d'une analyse de facies avec un algorithme de reconnaissance de forme supervise 
20 ou non. 



En reconnaissance des formes supervisee, il est necessaire de choisir des classes 
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d'apprentissage regroupant des evenements sismiques situes aux voisinages de puits dont 
l'interpr&ation g6ologique a 6te effectuee. Cette methodologie permet de classer 
T ensemble des evenements sismiques en fonction des caracteres g6ologiques observ6s aux 
puits. 

5 Dans notre cas, les connaissances g6ologiques sur les trois puits disponibles Wl, W2 

et W3, ainsi que les connaissances geologiques g6n6rales sur la zone d'etude nous amenent 
a considerer trois classes : une premiere classe correspondant aux puits Wl et W2 
repr^sentent des sables turbiditiques, une seconde classe correspondant au puits W3 
represente des boues turbiditiques et une troisieme classe repr6sente des h6mip61agites. Les 

10 puits Wl et W2 montrant des caract6ristiques geologiques proches, seuls les 6venements 
sismiques situ6s autour du puits W2 sont utilises, le puits Wl servant de test aveugle. La 
localisation du puits W3 6tant incertaine, les 6v6nements sismiques utilises pour 
caract6riser la classe 2 sont situ6s autour d'un puits virtuel d6fini & proximite de W3. 
Enfin, ne disposant d'aucun puits pour la classe 3, un puits virtuel est d6fini k partir de nos 

IS connaissances geologiques generates. La Fig.7 montre la localisation g6ographique des 
trois classes d'apprentissage. Chaque classe est constitute de 121 6v6nements sismiques 
(11x11). Ces evenements servent k calibrer une fonction de classement qui est ensuite 
appliqu6e k T ensemble des ev6nements sismiques. 

Les Fig.8A, 8B montrent respectivement les affectations obtenues et les affectations 
20 les plus probables, c'est-k-dire ayant une probability d' affectation sup6rieure ou 6gale k 
80%. Le puits Wl, utilis6 en test aveugle, est correctement attribu6 k la meme classe que le 
puits W2. Le puits W3 est localise en bordure des classes 1 et 2. La classe 2 se distribue le 
long d'une forme chenalisante. La carte de la figure 8B montre que les r£sultats sont 
stables : peu d'6v£nements presentent une mauvaise probability d' affectation. 

25 En reconnaissance des formes non supervisee, les classes d'apprentissage sont 

constitutes en utilisant les modes d'une fonction de density de probability multivariye 
calculee k partir de Tensemble des 6venements sismiques, dans un espace engendr6 par les 
huit premiers nouveaux attributs. Les evenements sismiques de plus forte density de 
probability constituant les pics de cette fonction de density sont syiectionnys pour ryaliser 

30 Tapprentissage de la fonction de classement. Cette fonction est ensuite appliquye k 
V ensemble des tenements sismiques considyrys pour les affecter a Tune des classes : les 
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Fig.9A, 9B montrent les affectations obtenues, ainsi que les affectations les plus probables, 
c'est-a-dire ayant une probability d'affectation sup&ieure ou 6gale k 80%. Les cartes 
obtenues sont tres proches des cartes des fig. 8A, 8B. La classe 3 de V analyse supervisee 
correspond aux classes 3, 4 et 5 de T analyse non supervis6e. Les puits Wl et W2 sont 
attribu6s & la meme classe et le puits W3 est toujours situ6 h. la frontiere entre les classes 3, 
et surtout 1 et 2. Cette convergence des resultats valide Interpretation de Tapproche 
supervisee. 
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REVENDICATIONS 

1) Mtthode pour compacter et filtrer des evenements geophysiques lisibles sur 
plusieurs jeux d'enregistrements gtophysiquesdisponibles, obtenus par exploration d'une 
zone souterraine, caracterisee en ce qu'elle comporte : 

5 - la r6partition de ces tvtnements geophysiques dans un certain nombre p de families de 
variables geophysiques ayant chacune une signification physique particuliere ; 

- la formation par combinaison des dites variables sismiques, de variables synthttiques 
en nombre tr&s inftrieur au nombre des dites variables, que Ton obtient par 
construction d'une base vectorielle orthogonale dans chacun des p ensembles d ! analyse 

10 constitutes par les donntes de chacune des dites families, d'oil dtcoule la formation 

d'une base vectorielle orthonormee decrivant les p ensembles d' analyse ; et 

- Tutilisation de cette base vectorielle orthonormee pour filtrer et decrire les dits 
evtnements geophysiques. 

2) Mtthode selon la revendication 1, caracterise en ce que la formation des variables 
15 synthttiques est obtenue en combinant une technique d' analyse en composantes principales 

ou des mtthodes dtrivtes de diagonalisation d'une norme et une technique de regression 
simple. 

3) Mtthode selon la revendication 1 et 2, caracttriste en ce que Ton soumet k une 
rotation les variables synthttiques pour amtliorer les correlations entre variables initiates et 

20 variables synthttiques. 

4) Mtthode selon la revendication 1 ou 2 ou 3, caracttriste en ce que Ton forme des 
variables synthttiques k partir de families constitutes de cubes de traces sismiques obtenus 
k Tissue d'une campagne sismique 3D, chacun correspondant k un meme deport ou offset. 

5) Mtthode selon la revendication 1 ou 2 ou 3, caracttriste en ce que Ton forme des 
25 variables synthttiques k partir de families constitutes de cubes de traces sismiques obtenus 

k Tissue d ! une campagne sismique 3D, chacun correspondant k un meme angle d' incidence. 

6) Mtthode selon la revendication 1 ou 2 ou 3, caracttriste en ce que Ton forme des 
variables synthttiques k partir de families constitutes de cubes d'attributs sismiques 
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obtenus par traitement des amplitudes de traces sismiques obtenues k Tissue d'une 
campagne sismique 3D. 

7) M6thode selon Tune des ievendications prec6dentes, caract6ris6e en ce que Ton 
foime des variables synthetiques k partir de families constituees de cubes de traces 

5 sismiques obtenus par des explorations sismiques successives de la zone. 

8) M6thode selon Tune des revendications pr6cedentes, caracterisee en ce que Ton 
forme des variables synthetiques k partir de families constituees de cubes d'attributs 
sismiques form6s par traitement des amplitudes de traces sismiques obtenues par des 
explorations sismiques successives de la zone. 

10 9) Methode selon la revendication 1 ou 2 ou 3, caract£ris6e en ce que Ton forme des 

variables synth6tiques k partir de families constitu6es de mesures d'attributs diagraphiques 
de nature physique diff6rente. 

10) Methode selon la revendication 9, caract6ris6e en ce que Ton forme des variables 
synthetiques k partir de families de mesures de resistivite, 

15 1 1) Methode selon la revendication 9, caract6risee en ce que Ton forme des variables 

synthetiques a partir de families de mesures li6es k la radioactivite des terrains. 

12) Methode selon la revendication 9, caracterisee en ce que Ton forme des variables 
synthetiques k partir de families de mesures li6es k la porosite des terrains. 
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